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一 种 多 目标 不 等 面积 设施 布局 问题 的 启发 式 算 法 
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摘 要 : 多 目标 不 等 面积 设施 布局 问题 《UA-FLP) 是 将 一 些 不 等 面积 设施 放置 在 车 间 内 进行 布局 ， 要 求 优 化 多 个 目标 
并 满足 一 定 的 限制 条 件 。 以 物料 搬运 成 本 最 小 和 非 物 流 关 系 强 度 最 大 ， 0 目标 优化 模型 ， 并 提出 一 种 
启发 式 算 法 进行 求解 。 算 法 采用 启发 式 布局 更 新 策略 更 新 构 型 ， 并 通过 结合 基于 自 适应 步 长 梯度 法 的 局 部 搜索 机 制 和 
启发 式 设施 变形 策略 来 处 理 设 施 之 间 的 干涉 性 约束 。 为 了 得 到 问题 的 Pareto ae 提出 了 基于 Pareto 优化 的 局 部 
搜索 和 基于 小 生境 技术 的 全 局 优化 方法 。 通 过 两 个 典型 算 例 对 算法 性 能 进行 测试 ， 实 验 结果 表明 ， 所 提出 的 启发 式 算 
法 是 求解 多 目标 UA-FLP 的 有 效 方法 。 
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Heuristic algorithm for unequal area facility layout problem with multiple objectives 


~ Liu Jingfa* ?, Liu Jun*? 
(a. Jiangsu Engineering Center of Network Monitoring, b. School of Computer & Software, Nanjing University of Information 
Science & Technology, Nanjing 210044, China) 


Abstract: The multi-objective unequal area facility layout problem (UA-FLP) is the problem of placing departments with 


different areas in a facility so that departments satisfy some given objectives and constraints. This paper built a multi-objective 
Oo optimization model based on minimizing the material handling costs and maximizing the closeness rating scores, and proposed 
: a heuristic algorithm (HA) to solve the multi-objective UA-FLP. This paper proposed the heuristic layout updating strategy to 
update the layouts, and used the combination of the local search based on the adaptive gradient method and the heuristic 
department deformation strategy to deal with the interface between any two different departments. The proposed heuristic 
algorithm combines the local search based on the Pareto optimization and the global search based on the niche technology to 
obtain Pareto-optimal solutions of the problem. Using two representative instances from the literature to test the performance 
of the algorithm, the experimental results show that the proposed heuristic algorithm is an effective method for solving the 
multi-objective UA-FLP. 


Key words: facility layout problem; heuristic algorithm; multi-objective optimization; Pareto optimal; niche technology 


了 一 些 求解 算法 。Komarudin 等 人 外 提出 了 基于 切片 树 特征 的 

蚁 群 系统 算法 ， 并 结合 了 几 种 局 部 搜索 方法 以 提高 算法 的 搜索 

设施 布局 问题 (facility layout problem,FLP) 是 指 将 若干 个 待 性 能 。Goncalve 等 人 9 建立 了 一 种 线性 规划 模型 ， 并 提出 了 一 

6 局 的 设施 放置 在 给 定 的 车 间 内 ， 在 满足 一 定 的 约束 条 件 下 ， 个 带 偏向 的 随机 密 钥 的 遗传 算法 。Liu 等 人 外 提出 了 一 种 基于 序 
优化 一 个 或 多 个 目标 函数 。 合 理 的 设施 布局 能 够 帮助 企业 提高 。 列 对 特征 的 模型 ， 并 通过 结合 遗传 算法 和 混合 整数 规划 方法 进 
生产 效率 ， 降 低 生 产 成 本 ,缩短 产品 周期 。 因 此 ， 车 间 设 施 布 。 行 求解 。Vitayasak 等 人 外 提出 了 一 种 基于 遗传 算法 的 回 溯 搜索 
局 问题 在 制造 业 中 具有 重要 的 研究 意义 。 算法 求解 UA-FLP， 其 中 回 济 搜 索 算 法 是 一 种 基于 群体 的 迭代 
近年 来 ， 一 些 学 者 对 设施 布局 问题 进行 了 广泛 的 研究 , 但 。 ”进化 算法 。Paes 等 人 中 提出 一 种 结合 分 解 策略 的 遗传 算法 ， 该 

是 大 多 数 文献 是 针对 相等 面积 设施 的 布局 问题 中 进行 讨论 , 在 。” ”算法 在 经 典 Benchmark 实例 上 进行 了 测试 ， 并 取得 了 高 质量 的 
实际 生产 应 用 中 ， 相 等 面积 的 设施 不 符合 实际 车 间 布 局 情况 。 布局 解决 方案 。 Asl 等 人 [提出 了 一 种 改进 的 粒子 群 算法 , 该 算 
前 ， 对 于 不 等 面积 设施 布局 问题 (UA-FLP)， 一 些 学 者 法 应 用 了 两 种 局 部 搜索 方法 和 设施 交换 方法 来 提高 布局 方案 的 
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质量 ， 防 止 局 部 最 优 解 。 印 胜 海 等 人 久 提 出 了 一 种 改进 的 遗传 
算法 ， 通 过 对 螺杆 制造 车 间 布 局 设计 方案 进行 优化 ， 实 验 结果 
证 明了 所 提出 的 算法 的 有 效 性 。 雇 源泉 等 人 外 提出 了 一 种 采用 
顺序 编码 方式 的 遗传 算法 ， 对 车 间 布 局 进行 优化 设计 以 减少 物 


流 搬运 费用 。 


过 


尽管 以 上 方法 对 求解 UA-FLP 取得 了 不 错 的 进展 ， 但 优化 四 1 设施 布局 示意 图 

的 目标 主要 和 集中 在 以 物料 搬运 成 本 为 主 的 单个 目标 优化 ， 忽 略 。” 1.2 数学 模型 
了 其 他 目标 的 实现 ， 比 如 空间 利用 率 、 非 物流 关系 强度 等 。 对 a) 目标 函数 
于 多 目标 UA-FLP， 一 般 使 用 线性 加 权 求 和 的 方法 将 多 目标 转 minf0O -党 3 Cpa, 而 
化 成 单 目 标 进行 求解 ， 但 是 线性 加 权 法 中 权重 系数 难以 合理 设 BI 
各 王 的 局 限 1 [10] 。 t 方 汗 3 中 站 N NWN 

A 对 此 ，Pareto 优化 方法 成 为 一 和 i te 而 
的 多 目标 求解 方法 。 徐 立 云 等 人 00 针 对 UA-FLP， 提 出 一 种 自 证 
适应 遗传 算法 ， 通 过 对 其 交叉 和 变异 算 子 进行 非 线性 化 处 理 ， 目标 函数 式 (1) 表 示 最 小 化 物料 搬运 费用 (MH Cost) ， 其 中 


克服 了 算法 过 早 收敛 和 陷入 局 部 最 优 的 缺点 。Ripon 等 人 09 提 ” ”Ci 为 设施 Di 与 Dj 之 间 的 单位 物料 单位 距离 搬运 费 ; 方 为 设施 
出 了 一 种 基于 进化 方法 的 自 适应 可 变 邻 域 搜索 的 算法 ， 以 最 小 。” Di 和 Dj 之 间 物 料 搬 运 频 率 ; di 为 设施 D; 和 Dj 之 间 的 物料 搬运 
物料 搬运 成 本 和 最 大 非 物流 关系 强度 为 优化 目标 ， 实 验 结果 表 ”距离 , 即 设施 中 心 点 之 间 的 欧 氏 距离 。 目 标 函 数 式 (2) 表示 最 大 
明 算法 具有 较 好 的 性 能 。 Kumar 等 人 03 提 出 一 种 基于 整数 线性 。” ”化 非 物流 关系 强度 (CR Score) ， 方 为 设施 之 间 的 非 物流 密切 度 
规划 和 变 邻 域 的 方法 求解 多 目标 UA-FLP。 最 近 ，Hunagund 等 。” 值 (如 表 1 所 示 〉。 如 果 设 施 之 间 存 在 共同 的 边界 ， 则 根据 设 
人 104 针对 动态 的 UA-FLP， 提 出 了 一 种 基于 弹性 区 带 结构 的 模 。 ” 施 间 的 密切 度 等 级 ， 计 算 密 切 度 值 ， 否 则 密切 度 值 为 0。 设 施 
拟 退 火 算法 。Liu 等 人 0 以 物料 搬运 费 、 总 邻接 值 和 车 间 利 用 ”之 间 非 物流 关系 的 值 越 大 ， 表 明 它 们 之 间 非 物流 关系 越 密切 。 


一 、 


率 为 优化 目标 ， 提 出 了 一 种 结合 目标 空间 分 割 方法 的 多 目标 粒 b) 约束 条 件 
子 群 化 算法 。 五 站 万 = 办 i et{l;2,., NBi¥ } (3) 
本 文 以 最 小 物料 搬运 费 和 最 大 非 物流 关系 强度 为 优化 目 x—0.5l >0, ie{l,2,...,N} (4) 
标 ， 通 过 应 用 基于 Pareto 优化 的 局 部 搜索 和 基于 小 生境 技术 的 X+0.5l SL, ie{l,2,...,N} (5) 
全 局 优化 搜索 ， 以 及 多 种 启发 式 策略 相 结 合 ， 为 解决 UA-FLP y—0.5w >0, ie{l,2,...,N} (0) 
提出 了 一 种 启发 式 算法 (HA) .通过 对 文献 中 的 两 个 代表 性 算 例 y+0.5w SW, ie{l,2,...,N} (7) 
进行 计算 ， 实 验 结果 表明 ， 所 提出 的 启发 式 算法 是 一 种 求解 多 maxfi mw < Re C8) 

标 UA-FLP 的 有 效 算 法 。 min{l, w} 

1 ”问题 描述 与 数学 模型 er Ee. 
其 中 : 和 表示 设施 D; 和 疡 的 内 部 空间 ，w 和 * 分 别 表 示 设 
1.1 问题 描述 施 Di 最 大 纵横 比 的 值 和 最 小 边 长 的 值 。 约 束 式 (3) 表 示 设 施 D; 


段 设 所 有 设施 及 车 间 的 形状 均 为 算 形 ， 设 施 的 长 宽 可 变 ， 和 Dj 之 间 互 不 嵌入 ; 约束 式 (4)~(7) 表 示 设 施 在 车 间 内 部 且 与 车 
但 其 面积 固定 ， 设 施 的 面积 之 和 等 于 或 小 于 车 间 面 积 。 多 目标 。 间 边 界 之 间 不 存在 嵌入 ; 约束 式 (8) 表 示 所 有 设施 满足 最 大 纵横 


UA-FLP 可 描述 为 : 给 定 矩 形 车 间 的 长 宽 ， 所 有 待 布设 施 的 面 比 限制 ;约束 式 (9) 表 示 所 有 设施 满足 最 小 边 长 限制 。 

积 ， 设 施 之 间 的 物流 量 ， 单 位 物料 搬运 费用 以 及 设施 之 间 的 非 表 1 设施 关系 密切 度 分 类 

沟 流 关系 强度 ， 要 求 给 出 布局 方案 使 得 布局 的 物料 搬运 费 最 小 密切 度 等 级 等 级 代号 密切 度 值 ry 

且 非 物流 关系 最 大 ， 并 满足 以 下 约束 条 件 ，a) 所 有 设施 必须 放 绝对 必要 A 6 

置 在 车 间 内 部 ， 设 施 之 间 不 能 存在 嵌入 ， 且 设施 与 车 间 边 界 之 特别 重要 E 5 

间 不 存在 嵌入 ; b) 所 有 设施 满足 最 小 边 长 或 者 最 大 纵横 比 限制 。 重要 I 4 
如 图 1 所 示 , 假设 笛 卡 尔 坐 标 系 的 原点 位 于 秆 形 车 间 据 的 一 般 3 

左下 角 , 工 和 W 分 别 是 矩形 车 间 邓 的 长 和 宽 ， Li 和 wi 分 别 表示 不 重要 2 

设施 Di 的 长 和 宽 。 设施 的 放置 方式 有 水 平 或 垂直 两 种 方式 。 如 不 能 靠近 1 

果 设 施 的 长 边 与 x 轴 平 行 ， 则 为 水 平 放 置 ， 如 果 设 施 的 长 边 与 


y 轴 平行 ， 则 为 垂直 放置 。 一 个 布局 〈 或 称 构 型 ) 用 
X=(ODyDIWDt2y2022 CNVYyNIwWwN 表 示 ， 其 中 (人 表示 设施 
Di 的 中 心 点 坐标 ，N 是 车 间 中 设施 的 数目 。 
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(a) 设施 与 车 间 边 界 在 x 轴 方 向 相交 


ee 
ci L = 一 一 一 jpr 


lir=|cipr 


(b) 设施 与 车 间 边 界 在 y 轴 方 向 相交 


Dp, 、 


lir=|dip#l 
(c) 设施 与 车 间 边 界 在 x 轴 和 >》 轴 方 向 同时 相交 
图 2 设施 Di 与 车 间 边 界 之 间 艇 入 示意 图 
借鉴 拟 物 思 想 09， 将 所 有 设施 内 部 和 车 间 外 部 均 想 象 为 光 
滑 的 弹性 实体 。 根 据 弹 性 力学 的 原理 ， 当 两 个 不 同 设施 之 间 存 
在 干涉 或 者 设施 与 车 间 边 界 之 间 存 在 干涉 时 ， 它 们 之 间 的 挤 压 


弹性 势能 正比 于 妃 入 深度 的 平方 527， 所 以 整个 系统 的 弹性 势能 
可 定义 为 
E(X)= Do + (10) 


其 中 : 4 是 弹性 系数 ， 本 文 设 置 /=1; E(X) 表 示 布 局 的 总 的 弹性 
势能 ;i 表示 设施 Di 和 疡 之 间 的 嵌入 深度 ， 兢 表示 设施 Di 与 
车 间 边 界 之 间 的 嵌入 深度 。 如 果 设 施 Di 与 Dj 互相 嵌入 ， 则 符 
入 深度 定义 为 将 设施 D; 和 Dj 变 为 不 相交 状态 所 需要 在 x 轴 或 
者 y 轴 方向 上 移动 的 最 短 距离 071。 如 果 设 施 Di 与 车 间 边 界 相 交 ， 
可 分 为 三 种 情况 a) 设施 Di 与 车 间 边界 在 x 轴 方 向 相交 (如 图 
2(a) 所 示 ) ; b) 设 施 Di 与 车 间 边 界 在 y 轴 方 向 相交 (如 图 2(b) 
所 示 ) ; 9 设施 Di 与 车 间 边界 在 x 轴 方 向 与 y 轴 方 向 同时 相交 

(如 图 2(c) 所 示 〉 ， 则 杉 入 深度 分 别 为 lir=|bip#l, |cips| 和 |dip#l， 
其 中 |AB| 表 示 顶 点 A 与 顶点 B 之 间 的 欧 氏 距离 。 


2 ”算法 设计 

对 于 多 目标 UA-FLP， 为 了 得 到 一 组 Pareto 最 优 解 ， 本 文 
提出 一 种 启发 式 算法 (HA)。 首 先 随 机 产生 一 组 规模 为 n 的 初始 
布局 ( 构 型 ), 针对 每 种 布局 , 提出 一 种 启发 式 布局 更 新 策略 来 
更 新 当前 布局 ， 再 运用 结合 了 基于 自 适 应 步 长 的 梯度 法 和 启发 
式 设施 变形 策略 的 布局 合法 化 机 制 , 使 每 个 布局 满足 约束 条 件 ; 
然后 提出 基于 Pareto 优化 的 局 部 搜索 和 基于 小 生境 技术 的 全 局 


Chin: 
面积 


刘 景 发 ， 等 : 一 种 多 目标 不 等 


| 内 


弹性 势能 CP; 最 大 的 设施 应 进 行 重新 放置 ， 其 中 CP; 的 计算 方 
式 见 公式 (11)， 其 中 hj; 是 设施 Di 的 面积 。 

N 

pa 有 二 


CP= 刁 二 
1 A 3» 


i NN (IL) 


如 果 当 前 布局 是 合法 布局 ， 则 挑选 出 相对 物料 搬运 费 MC 


MC = ,11,2,.,N. (12) 


为 了 增加 布局 的 多 样 性 ， 设 施放 置 的 方式 为 :将 被 挑选 设 
施 的 中 心 随 机 放置 在 车 间 空 白 区 域内 2N 次 (如 果 车 间 内 没有 
空白 区 域 ， 则 将 被 挑选 设施 的 中 心 随 机 放置 在 车 间 内 2N 次 )， 
且 每 次 放置 均 包 含水 平 放 置 和 垂直 放置 两 种 方式 ， 这 样 一 共 得 
到 4N 个 新 构 型 。 
2.2 布局 合法 化 
当 通 过 启发 式 布局 更 新 策略 产生 4N 个 新 构 型 后 ， 这 些 构 
型 不 一 定 满足 约束 条 件 ， 构 型 中 可 能 存在 设施 之 间或 设施 与 车 
间 边 界 之 间 相 互 嵌 入 。 本 文通 过 交替 执行 局 部 搜索 机 制 和 启发 
式 设施 变形 策略 对 不 合法 的 构 型 进行 布局 合法 化 操作 ， 使 得 所 
有 构 型 满足 约束 条 件 。 
2.2.1 局 部 搜索 机 制 
当 通 过 启发 式 布局 更 新 策略 得 到 一 组 新 布局 后 ， 为 了 得 到 
布局 附近 弹性 势能 更 低 的 布局 ， 本 文 引 入 一 种 基于 自 适 应 步 长 
梯度 法 (GMDQ8 的 局 部 搜索 机 制 , 使 设施 之 间或 设施 与 车 间 边 界 
之 间 的 嵌入 部 分 由 于 弹性 力 的 作用 ， 相 互 推 开 ， 从 而 获得 合法 
布局 。 
GM 又 被 称 为 最 速 下 降 法 ， 其 搜索 方向 为 函数 值 下 降 最 快 


的 负 梯 度 方向 。 在 GM 和 迭代 过 程 中 ， 对 于 一 个 构 型 六， 先 计算 
X! 在 x 轴 方 向 和 y 轴 方 向 的 梯度 ， 从 而 得 到 新 构 型 Z2= XI -hh x 


VCECXD))， 其 中 大 为 迭代 步 长。 如果 新 产生 构 型 夭 的 能 量 E(X”) 
比 X! 的 能 量 ECX') 高 ， 则 表明 步 长 h 太 大 ， 通 过 将 h 乘 以 系数 
0.8 以 减少 步 长 ， 否 则 维持 步 长 不 变 。 为 了 缩短 GM 的 运行 时 
间 ， 一 种 “提前 逃离 ”策略 被 应 用 于 GM， 即 如 果 新 产生 构 型 
XX 与 构 型 如 的 能 量 很 接近 (|E(X”) - ECXCDJ|<10) 且 新 构 型 作 的 能 
量 EC 大 于 1， 则 认为 构 型 天 不是 一 个 有 “前 景 ”的 构 型 ， 
GM 提前 结束 ， 这 样 可 以 避免 由 于 GM 步 长 过 小 导致 搜索 时 间 
过 长 的 问题 。 如 果 新 产生 的 构 型 大 的 能 量 EC2) 小 于 1， 则 继 


搜索 的 两 种 搜索 方法 寻找 最 优 布局 ， 最 后 给 出 所 提出 启发 式 算 

法 的 具体 步 又 。 

2.1 启发 式 布局 更 新 策略 

在 提出 的 HA 中 ， 每 次 迭代 都 需要 更 新 布局 ， 根 据 车 间 布 

局 的 特点 ， 本 文 提出 一 种 启发 式 布 局 更 新 策略 ， 对 设施 进行 挑 

选 和 放置 。 
如 果 当 前 布局 是 非法 的 布局 ， 则 从 布局 中 挑选 出 相对 挤 压 


续 执 行 梯度 法 , 直到 构 型 的 能 量 满足 E(X?)< 8，, 或 者 迭代 步 长 有 h 
小 于 最 小 迭代 步 长 hwn， 其 中 s 是 一 个 合法 构 型 必须 达到 的 弹 
性 势能 的 最 小 值 。 事 实 上 ，h，hin 和 均 是 通过 反复 实验 得 到 
的 经 验 值 ， 当 he[0.1,10]，hnin EE[10”5, 10-3], EE[10-35, 10-5] 时 ， 
本 文 启发 式 算法 得 到 的 计算 结果 没有 什么 差别 ， 且 均 能 得 到 一 
个 合法 构 型 。 因 此 ， 本 文 设置 h=1，hwin=104，&e=10-20。 
适应 步 长 的 梯度 法 GM 的 详细 步骤 如 下 : 
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a) 设 


参数 有 = 1, hmin = 10"4, s=10-20 ; 


b) 计 算 当 前 构 型 X! 的 能 量 函 数 EX1) 的 x 轴 和 y 轴 方 向 的 


梯度 : 


VE(CXD = 


,一 一 ,一 一 
Bo" Or 


得 到 新 的 构 型 X= X! -hx VCECXD); 


c) 如 果 ECX2) >1， 则 


( 


a) 如 果 E(X?) > E(X!)， 则 令 hh=hx0.8; 


(b) X=X2, =X -hxV EC): 


( 


d) 如 果 


( 


c) 如 果 |E(X?”) - EC > 104， 则 转 步 又 (a); 
EX”) <1, 则 


a) 如 果 ECXZ2) > E(X!)， 则 令 及 = 有 hx 0.8; 


(b) X=X2, X=X1-hxV ECX): 


(c) 如 果 E(X?)> 8 且 有 > jnin， 则 转 到 步骤 (a); 
e) 输 出 和 大. 
2.2.2 启发 式 设施 变形 策略 


到 平衡 ， 
法 化 。 


在 执行 GM 过 程 中 ， 有 可 能 产生 “ 卡 壳 ”现象 ， 即 虽然 构 
型 中 存在 设施 嵌入 ， 但 是 嵌入 设施 的 受 力 在 两 个 不 同方 向 上 达 
此 时 无 法 通过 执行 GM 将 嵌入 的 设施 推 开 来 使 布 
于 设施 的 面积 是 固定 的 ， 在 满足 约束 条 件 式 (8) 和 (9) 
的 前 提 下 ， 适 当 改 变 和 矩形 设 施 的 长 宽 ， 可 以 得 到 合法 布局 。 因 
此 在 执行 完 GM 后 ， 提 出 一 种 启发 式 设施 变 


局 合 


策略 ， 来 使 布局 


启发 式 设 施 变 形 策略 过 程 如 下 : 首先 挑选 出 所 有 与 车 间 边 
界 嵌 入 的 设施 ， 按 照 设 施 与 车 间 边 界 嵌 入 深度 
进 


小 到 大 的 顺序 
操作 ， 即 将 设施 与 车 间 边 界 嵌 入 的 部 分 移动 到 设 


设施 变形 
抽空 白 的 


x 域内 ,同时 保证 变形 后 的 设施 仍然 是 矩形 设施 ， 


且 满 足 约束 式 (8) 和 (9) ## 
音 按照 嵌入 深度 从 小 到 大 的 顺序 依次 挑 
形 操作 ， 将 嵌入 部 分 转移 到 设施 周转 


次 ， 对 于 设施 之 间 相 互 嵌 入 的 情况 ， 
选 设施 进行 类 似 的 变 
的 区 域 。 判 断 


一 个 设 


施 是 否 适合 进行 变形 操作 的 依据 是 : 
局 系统 总 的 弹性 


设施 的 变形 ， 


人 下 

在 执行 完 设施 变形 后 ， 布 
势能 是 否 下 降 ? 如 果 弹 性 势能 下 降 ， 则 进行 该 
否则 不 执行 此 设施 的 变形 。 如 果 设 施 周 围 没 有 足 


够 的 空 


区 域 容纳 嵌入 的 部 分 ， 则 先 移动 嵌入 的 一 部 分 到 笃 
区 域 ， 然 后 将 剩余 的 藤 入 部 分 与 其 他 设施 与 其 嵌入 部 分 之 和 移 


动 到 此 设施 周围 其 他 的 空白 区 域内 。 如 果 执 行 完 启发 式 设施 变 
形 策略 以 后 ， 车 间 中 仍然 存在 设施 之 间或 设施 与 车 间 边 界 之 间 


的 嵌入 ， 或 设施 周围 存在 邱 


阶 ， 则 再 


次 执行 GM， 使 得 设施 之 


间 更 加 紧凑 ， 同 时 使 得 布局 能 成 为 满足 约束 的 合法 布局 。 
下 面 给 出 一 个 典型 的 例子 来 闭 述 设施 变形 的 过 程 ( 如 图 3 


所 示 )。 为 了 提高 阅读 性 , 将 车 间 空 


x 域 用 阴影 表示 。 假设 图 


3(a) 是 执行 完 GM 之 后 的 布局 状态 ， 其 中 Ds 


的 边界 ，D21 


嵌入 De 和 和 车间 
联 入 D1 和 D3。 按照 启发 式 设施 变 形 策略 的 原则 ， 


与 D6 的 嵌入 部 分 转移 到 Ds 下 方 包含 气 
(图 3(d)); 最 后 再 移动 D2 与 Di 
或 的 矩形 区 域 (图 3(e)), 在 此 过 程 中 布 
变 小 ， 但 执行 


A 


X 
不 断 


2.3 基于 Pareto 
对 多 
初始 构 型 ， 然 后 对 于 每 个 构 型 Xi(ie{1,2,.…, n})， 


式 布局 更 新 策略 和 布局 合法 化 操作 ， 得 到 一 组 由 


在 旬 


首先 挑选 与 车 间 边 界 嵌入 的 设施 Ds 进行 变 
间 赎 入 的 部 分 转移 到 Ds 上 方 
3())。 然 后 挑选 设施 之 间 相 互 嵌 入 的 设施 ， 按 照 嵌 入 深度 
到 大 的 顺 
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刘 景 发 ， 等 : 一 种 多 目标 不 等 面积 设 问题 的 启发 式 算 ; 


多 操作 ， 将 Ds 与 车 


全 
包 售 至 


又 域 的 矩 


x 域内 (图 


Rs 


又 域 的 矩形 


3(c) ); 


x 域内 (图 
由 小 


序 对 设施 进行 变形 操作 。 先 将 D2 与 D3 的 嵌入 部 分 转 
移 到 D? 左 侧 包 含 空 


其 次 将 Ds 


的 散 入 部 分 到 


完 所 有 设施 依次 变 


操作 后 ， 


区 域 的 矩形 区 


忆 右 侧 包 含 空 


局 整体 的 弹性 势能 


D5 与 D7 仍然 


嵌入 ， 再 次 执行 梯度 法 ， 最 终 得 到 合法 布局 ， 如 图 (3G) 
D; D; 
YG D. Do 
~ ~ 
D D 
器 
ZZ 人 ZZ 
六 Di 
D; D; 
Ds Ds 
2 9 
(a) 初始 布局 (b) Ds 变形 
Ds 六 
Ds Ds 
D; D; 
D D 
Ds Ds 
D 区 D 
ji D, 
D; D; 
D D 
2 2 
(c) D; 变形 (d) Ds 再 次 变形 
D; 家 
Do D 
D D 
D D 
D5 Ds 
D D 
Ds Ds 
D 
Ds Ds Db 
MM 
(e) D2 再 次 变形 (f) 再 次 执行 GM 
图 3 ”一 种 典型 的 启发 式 设施 变形 策略 执行 过 程 示意 图 


支配 的 局 部 搜索 


标 UA-FLP 的 HA 算法 中 ， 首 先 随机 产生 n 个 


通过 执行 启发 


4N 个 布局 组 


成 的 试验 性 构 型 库 , 记 为 TC。 在 UA-FLP 中 , 一 个 解 对 应 一 个 
布局 或 构 型 。 通 过 比较 当前 解 入 与 其 对 应 的 试验 性 构 型 库 TC 


中 解 吕 =1, 2,.…., 4N) 之 间 的 Pareto 支配 关系 00 来 决定 TC 中 
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解 的 质量 函数 的 值 以 及 当前 解 的 搜索 方向 。 如 果 试 验 性 构 型 
中 的 解 咏 是 不 可 行 的 解 ， 则 说 明 解 态 对 当前 解 总 的 搜索 方 
没有 帮助 ; 如 果 解 名 是 可 行 解 且 支 配 XX, 则 说 明 解 为 有 利于 当 
前 解 X 朝 着 Pareto 前 沿 的 方向 迭代 , 因此 解 的 质量 函数 的 值 也 
就 越 大 。 在 迭代 1 步 时 ， 当 前 解 X 相 对 于 TC 中 叉 的 质量 函数 
6(D 的 定义 如 下 : 


Pei 


人，X ,为 非 可 行 解 ; 
Po 1， 人 已 支 配 X ij; 
”| 42，X 为 可 行 解 且 X 与 X ,为 非 支配 关系 ; 
如，XX; 为 可 行 解 且 关 ,支配 X,. 
其 中 : je{1, 2…, 4N}。 参数 4, ?2,13,44 是 质量 函数 的 4 个 等 级 
且 机 < 和 p<h3<h4。 考 虑 到 寻 优 方向 除了 与 解 的 质量 有 关 之 外 ,也 
与 解 之 间 的 距离 有 关 。 实 验 性 构 型 库 中 解 的 质量 函数 值 越 大 ， 
与 当前 解 的 距离 越 近 ， 则 被 选 为 当前 解 寻 优 方向 的 概率 越 大 ， 
故 对 于 TC 中 的 每 个 实验 性 构 型 入， 选择 的 概率 定义 如 下 : 
0,(1)5,(7) 
D0,05,0) 
其 中 : 8(D) 是 算法 迭代 1t 步 时 大 的 质量 函数 ,6(D=1/f4di, 其 中 fds 
为 当前 解 成 与 试验 性 构 型 库 中 总 之 间 的 距离 ， 计 算 公 式 如 下 : 


(13) 


P.O) 沁 1,2,...,4N (14) 


fa = GAD) AX) (15) 


其 中 : je{1, 2,..…., 4N}，m 是 目标 函数 的 个 数 。 
2.4 基于 小 生境 技术 的 全 局 搜索 


仅仅 依靠 于 基于 Pareto 优化 的 局 部 搜索 方法 ， 算 法 可 能 陷 
入 局 部 极 小 值 ， 且 对 解 的 多 样 性 不 易 保 存 。 为 了 提高 算法 的 全 


局 搜索 能 力 ， 本 文 提 出 一 种 基于 小 生境 技术 的 全 局 搜索 方式 。 
设立 一 个 用 来 保存 当前 算法 找到 的 所 有 可 能 的 非 支配 解 的 外 部 
着 合 (文档 )， 记 为 BP。 当 算法 采用 基于 小 生境 技术 的 全 局 搜 
索 进行 寻 优 时 ， 挑 选集 合 BP 中 散布 最 为 稀疏 的 非 支配 解 Xi， 
作为 当前 解闷 的 搜索 方向 。 

假设 当前 BP 中 有 p 个 非 支配 解 Xi1, X2, … Xp， 对 于 当前 解 
Xi, 集合 BP 中 小 生境 数 最 小 的 非 支 配 解 敌 作为 XX 的 寻 优 方向 ， 
小 生境 数 niche( 有 ) 的 计算 公式 如 下 : 


niche(]) = by S(fd,), 1=1,2,.,7. (16) 


j=1,j#! 


其 中 : SV) 为 共享 函数 的 值 ， 其 计算 公式 如 下 : 


辐 fdj low, iffd; < owe; 
S(fd;)= (17) 
0, else. 
其 中 : Gynare 是 小 生境 半径 ，fiy 为 BP 中 两 个 解 尺 和 已 之 间 的 
距离 ， 其 计算 公式 参见 式 (15)。 
2.5 启发 式 算法 步骤 描述 

HA 的 具体 迭代 步骤 如 下 : 

a) 随 机 产生 n 个 初始 布局 。 设 Po= 0.6, t= 1， 最 大 友人 代步 
数 为 7。 

b) 对 于 每 一 个 布局 , 执行 GM 和 启发 式 设施 变形 策略 进行 
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由 和 


合法 化 操作 ， 得 到 一 组 尽 可 能 合法 的 布局 ， 称 为 构 型 库 ， 记 为 
C。 


c) 确 定 C 中 所 有 非 支配 解 ， 记 为 BP。 

d) 令 二 1。 

eo) 随机 产生 一 个 [0, 1] 内 的 随机 数 P, 当 P< Po 时 , 则 对 C 
中 的 当前 构 型 XX 采用 基于 小 生境 技术 的 全 局 搜索 ， 从 BP 中 挑 
选 出 最 为 稀 琉 的 非 支 配 解 , 记 为 X， 否 则 转 步 又 下。 

f) 对 于 C 中 的 当前 构 型 X%， 采 用 启发 式 布局 更 新 策略 ， 产 
E 4N 个 新 构 型 ， 对 每 个 新 构 型 执行 GM 和 启发 式 设施 变形 策 
略 ， 得 到 4N 个 尽 可 能 合法 的 构 型 ， 称 为 试验 性 构 型 库 ， 记 为 
TC。 采 用 基于 Pareto 优化 的 局 部 搜索 策略 ， 从 TC 中 挑选 出 一 
个 最 优 构 型 ， 记 为 总 。 

8g) 更 新 构 型 库 C， 即 令 X=X。 

h) 更 新 BP。 如 果 挑 选 的 可 行 解 X/ 对 于 BP 来 说 是 非 支配 
的 ， 则 将 构 型 入 加 入 BP， 并 删除 BP 中 被 恕 所 支配 的 解 ; 否 
则 BP 保持 不 变 。 

i) 令 i=i+1。 如果 i<n, 则 转 到 步 又 e); 否则 转 到 步骤 j)。 

j) 令 1= 1。 如 果 +t< 7T， 则 转 到 步骤 c); 否则 输出 BP， 
算法 结束 。 


3 ”仿真 实验 与 结果 分 析 
3.1 测试 算 例 


采用 文献 中 两 种 不 同 规模 的 普遍 使 用 算 例 对 所 提出 的 HA 
算法 进行 测试 。 算 例 08 由 Meller 等 09 首 次 提出 ， 将 8 个 设施 
放置 在 大 小 为 11.31x13 的 车 间 内 。 本 文 设 置 初始 构 型 的 个 数 n 
为 200， 人 迭代 步 数 了 为 100, 质量 函数 的 4 个 等 级 分 别 设 为 1= 
0.01，4=0.1，43=2，44=5; 算 例 SC30 由 Liu 和 Meller 册 提出 ， 
将 30 个 设施 放置 在 大 小 为 50x50 的 车 间 内 ， 本 文 设置 初始 构 
型 的 个 数 n 为 1000， 达 代 步 数 了 为 9000， 质量 函数 的 4 个 等 级 
分 别 设 为 1= 0.01，12=0.2，13=8，h=32。 对 于 每 个 算 例 ， 单 位 
距离 单位 物流 搬运 费 都 设 为 1。 关 于 物流 量 ， 设 施 的 面积 ， 密 
切 度 等 级 ， 最 大 纵横 比 约束 和 最 小 边 长 约束 等 详细 参数 设置 ， 
均 可 在 文献 外 中 找到 。 
3.2 ”实验 结果 与 分 析 

在 一 台 配 置 为 Intel Core 2 Duo, 2.94 GHz CPU 和 2.0 GB 
RAM 的 PC 上 , 利用 Java 语言 对 HA 算法 进行 运算 。 对 于 每 个 
算 例 , 分 别 独立 运行 HA 算法 30 次 , 对 HA 算法 得 到 的 目标 最 
好 值 (Best) 与 平均 值 (Avg) 这 两 个 方面 进行 评估 ,并 与 Ripon 等 [3 
提出 的 自 适 应 变 邻 域 搜索 算法 (VNS) 和 Ripon 等 2 提出 的 带 局 
部 搜索 的 进化 算法 (EA with LS) 和 不 带 局 部 搜索 的 进化 的 算法 
(EA without LS) 进 行 结果 对 比 ， 并 将 实验 结果 列 在 表 2 中 。 

根据 表 2 中 的 结果 可 知 ，HA 对 于 物流 搬运 费 (MH cost) 和 
非 物流 关系 强度 (CR score) 在 最 好 值 和 平均 值 上 都 优 于 或 等 于 
其 他 三 种 算法 。 对 于 算 例 08， 本 文 算法 相 比 于 三 种 对 比 算 法 中 
结果 最 好 的 VNS 方法 , HA 得 到 的 关于 目标 MH Cost 的 结果 改 
进 了 (202.72-202.29)/202.71=0.21%， 而 CR Score 的 值 相同 。 对 
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于 算 例 SC30， 本 文 HA 算法 的 结果 


上 VNS 算法 得 到 的 关于 目 


标 MH cost 的 结果 改进 了 (3707-3553.53)/3707=4.14%， 而 CR 


score 改进 了 (392-380)/392= 3.06%。 图 
得 到 的 算 例 08 和 SC30 


别 表 示 


算法 HA 


4 和 图 5 用 三 角形 点 分 
的 Pareto 最 优 解 集 。 


从 图 中 可 以 看 出 ， 
有 良好 的 多 样 性 ， 
表 2 四 种 不 


通过 HA 算法 得 到 的 Pareto 解 分 布 均匀 ， 


能 形成 明显 的 Pareto 前 端 。 


同 算法 得 到 算 例 08 和 SC30 的 结果 比较 


物料 搬运 费 非 物 流 关 系 强度 
(MH Cost) (CR Score) 
算 例 对 比 算法 
最 优 值 平均 值 最 优 值 平均 值 
(Best) (Avg) (Best) (Avg) 
O8 HA( 本 文 方法 ) 202.29 251.157 70 51.428 
VNSI03 202.72 270.125 70 42.328 
EA with LSP9 202.72 279.69 68 38.00 
EA without LSP9 202.72 292.602 64 32.267 
SC30 HA (本 文 方法 ) 3553.53 4101.84 392 287.48 
VNSI03 3707.00 4146.27 380 264.62 
EA with LSP9 3716.44 4186.04 370 249.834 
EA without LSP9 3740.75 4444.037 349 241.872 


为 了 进 


步 分 析 算 法 的 性 能 ， 本 文 分 别 给 出 了 两 个 算 例 目 


标的 最 优 值 (best) 和 平均 值 (avg) 随 着 迭代 步 数 增加 的 变化 情况 。 


的 兴 代 图 


图 6(a) 和 图 6(b) 分 别 是 算 例 O08 的 两 个 有 
， 图 7(a) 和 
cost 和 CR score 的 迭代 图 , 图 中 


标 MH Cost 和 CR score 
图 7(b) 分 别 是 算 例 SC30 的 两 个 目标 MH 
标 值 取 自 于 HA 算法 30 次 运 


行 的 结果 。 从 图 


波动 比较 大 , 但 是 随 着 迭 


定 ， 说 明了 算法 具有 较 好 的 收敛 速度 和 稳定 折 


中 可 以 看 出 ， 


两 组 算 列 的 目标 在 算法 迭代 初期 ， 
代步 数 的 增加 ，Best 和 Avg 都 趋 于 稳 
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(b) CR score 迭代 图 
6 算 例 08 目标 运 代 图 


MH Cost (minimize) 


Generations 


(a) MH cost 迭代 图 
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4 ”结束 语 
本 文 研究 多 目标 UA-FLP， 通 过 结合 GM 和 启发 式 设施 变 


形 策略 ， 提 出 了 一 种 


高 效 的 解决 设施 之 间 互 不 能 入 的 约束 处 理 


机 制 。 该 约束 处 理 方法 是 一 种 解决 〈 二 维 ) 车 间 设 施 布局 问题 


中 设施 之 间 不 租 入 的 


直接 简单 的 推 
过 将 基于 Pareto 的 局 


问题 ， 但 从 设施 变形 
施 布局 ， 对 于 圆 形 、 三 角形 、 以 及 不 规则 形状 的 任意 设施 变形 ， 
k 有 困难 ， 这 将 是 以 后 看 


有 效 方法 ， 同 时 也 能 推广 到 三 维 矩形 布局 


过 程 不 难看 出 : 该 方法 主要 适用 于 矩形 设 


究 的 方向 。 
部 搜索 和 基于 小 生境 的 全 局 搜索 两 种 搜索 


方式 相 结 合 ， 避 和 免 ] 


能 够 避免 算法 陷入 局 


通常 使 用 单一 搜索 方式 的 局 限 性 ， 这 不 仅 
部 最 优 ， 增 强 算法 的 全 局 搜索 能 力 ， 而 且 


录用 稿 


能 够 提高 算法 获得 Pareto 最 优 解 的 多 样 性 。 通 过 两 个 典型 的 自 
例 对 算法 性 能 进行 测试 和 分 析 ， 并 与 目前 文献 中 优秀 的 算法 进 
行 对 比 。 实 验 结果 表明 ， 本 文 提 出 的 HA 算法 对 于 解决 多 目标 
A-FLP 是 一 种 有 效 的 方法 ， 在 解 的 收敛 性 和 人 解 的 多 样 性 等 方 
面 均 具有 明显 的 优势 ， 能 够 为 企业 提供 有 效 的 布局 解决 方案 ; 
同时 ， 所 提出 的 HA 算法 也 能 推广 应 用 于 解决 其 他 多 目标 优化 


问题 。 
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